SACLAY-T10/062 
CERN-PH-TH/2010-097 



O 

(N 



(N 
(N 

in: 

(N 
> 

OS 
t> 
t> 

O 
O 



X 



Can multistate dark matter annihilation explain 
the high-energy cosmic ray lepton anomalies? 

Marco Cirelli 1 ' 2 , James M. Cline 1 - 3 
1 CERN Theory Division, CERN, Case CO 1600, CH-1211 Geneve, Switzerland 

2 Insitut de Physique Theonque, CNRS URA 2306 & CEA/Saclay, F-91191 Gif-sur-Yvette, France 

and 
3 Physics Department, McGill University, 3600 University Street, Montreal, Quebec, Canada H3A 2T8 
e-mail: marco.cirelli@cea.fr, jcline@physics.mcgill.ca 
(Dated: 17 May 2010) 

Multistate dark matter (DM) models with small mass splittings and couplings to light hidden 
sector bosons have been proposed as an explanation for the PAMELA/Fermi/H.E.S.S. high-energy 
lepton excesses. We investigate this proposal over a wide range of DM density profiles, in the 
framework of concrete models with doublet or triplet dark matter and a hidden SU(2) gauge sector 
that mixes with standard model hypercharge. The gauge coupling is bounded from below by the DM 
relic density, and the Sommerfeld enhancement factor is explicitly computable for given values of the 
DM and gauge boson masses M, fi and the (largest) dark matter mass splitting 5Mi2- Sommerfeld 
enhancement is stronger at the galactic center than near the Sun because of the radial dependence 
of the DM velocity profile, which strengthens the inverse Compton (IC) gamma ray constraints 
relative to usual assumptions. We find that the PAMELA/Fermi/H.E.S.S. lepton excesses are 
marginally compatible with the model predictions, and with CMB and Fermi gamma ray constraints, 
for M = 800 GeV, p £ 200 MeV, and a dark matter profile with noncuspy Einasto parameters 
a iZj 0.20, r s ~ 30 kpc. We also find that the annihilating DM must provide only a subdominant 
( ,$, 0.4) component of the total DM mass density, since otherwise the boost factor due to Sommerfeld 
enhancement is too large. 

PACS numbers: 95.35.+d, 98.70.Sa, 12.60Cn 



I. INTRODUCTION 

Dark matter (DM) annihilating in the galaxy has 
been proposed as the explanation for a number of 
current experimental anomalies. The Payload for 
Antimatter-Matter Exploration and Light-nuclei Astro- 
physics (PAMELA) experiment has observed an excess in 
the positron fraction in the energy range 10—100 GeV [l|, 
and the Fermi Large Area Telescope (LAT) [2| observes 
an excess in e + + e - at energies up to 1 TeV, compared 
to the predictions of conventional diffusive models. The 
High Energy Stereoscopic System (H.E.S.S.) [3], too, ob- 
serves a steepening in the spectrum around a few TeV 
which is in agreement with the Fermi-LAT excess. On 
the other hand, PAMELA does not observe excess an- 
tiprotons. Therefore TeV-mass particles that can anni- 
hilate preferentially into leptons are the preferred DM 
candidates for explaining these anomalies 0t3|- 

Theories of leptophilic annihilating DM have come un- 
der increasing pressure from a number of complementary 
constraints, due to the radio and gamma ray (GR) emis- 
sions that should be produced by the annihilation pro- 
cess itself ('prompt' gamma rays) [9|-Tl2i| or by the high 
energy e , via synchrotron radiation in the galactic mag- 
netic field or via inverse Compton scattering on cosmic 
microwave background (CMB)_photons and the galactic 
radiation field (starlight) [13H20J - Such radiation could 
also partially reionize the universe after recombination, 
leading to distortions in the CMB [2lH25J|. Finally, rel- 
evant constraints are imposed by extragalactic gamma 



rays, produced by the annihilations in forming DM ha- 
los J22, IM El- Many models of annihilating DM are 
seemingly ruled out by overproducing such radiation. In 
particular, inverse Compton gamma rays produced in the 
galaxy would contribute to the GR spectrum observed by 
Fermi-LAT. 

A precise assesment of the strength of Fermi inverse 
Compton (IC) constraints depends however upon the 
details of the assumed DM density profile p(r). The 
strongest constraints come from close to the Galactic 
Center (GC); therefore they can be satisfied if p(r) is 
not strongly peaked as r — > 0, where r is the distance 
from the GC. Usually such dependence is illustrated by 
choosing a few different paramctrizations of the profile 
with varying cusp iness (e.g., the peaked Navarro-Frcnk- 
White profile [28[ as opposed to the cored isothermal [29| 
or Burkert profiles 30]). iV-body simulations and obser- 
vations of our own galaxy suggest a reasonable range of 
shapes for p(r), and it would be interesting to explore 
the dependence of the constraints on the characteristics 
of p(r) within such a range. One purpose of the present 
work is to fill this gap by quantifying the strength of 
the IC constraint for a wide range of DM density pro- 
files. We limit ourselves to the Einasto parametrization 
[3ll |32| , which we find to be useful for this purpose since 
it allows for continuous variation in the cuspiness and it 
has been used to fit the results of high- resolution A^-body 
simulations. 

The constraints also depend upon which leptons, and 
how many of them, are assumed to be produced in the 



DM annihilations. For example, models that produce 
only e + e _ are more strongly constrained, for a given DM 
mass M, than those that produce two e^ pairs, since 
the spectrum of the former is harder. This is a relevant 
distinction since some of the models that are most nat- 
ural from the particle physics point of view necessarily 
produce the 4e (or at least 2e plus an invisible hidden 
sector gauge boson) rather than the 2e final state. Like- 
wise models that produce muons are more strongly con- 
strained than those that produce only electrons. More- 
over models using the gauge kinetic mixing portal pro- 
duce an admixture of e , fj, and it rather than just a 
single kind of final state. Failing to recognize these dis- 
tinctions, one could be misled into thinking that the gen- 
eral class of models is excluded, even though there may 
exist specific examples that satisfy all the constraints. 

Thus another goal of the present work is to explore 
in detail the viability of a very specific set of annihilat- 
ing DM models, which are also theoretically motivated 
[33|-|36|. Namely we consider DM %i transforming under 
a nonabelian hidden sector gauge symmetry, taking the 
simplest example of SU(2) and the lowest representations 
(doublet or triplet) for the DM. One of the dark gauge 
bosons Bi is assumed to mix with standard model hy- 
percharge (and thus the photon) through a dimension-5 
gauge kinetic mixing operator. This model has a num- 
ber of appealing features. The DM annihilates into dark 
gauge bosons, \X ~^ BB. By assuming the gauge sym- 
metry breaks below the GeV scale, so the Bi's have mass 
(j, ;$ 1 GeV, it is assured that only light leptons or n ± are 
produced by the gauge boson decays, and no antiprotrons 
that would be in conflict with the PAMELA observations. 
The annihilation cross section is naturally enhanced by 
Sommerfeld effect [33, |37| in this model, allowing it to 
have the correct relic density in the Early Universe, even 
though the cross section in the galaxy must be larger 
by a factor of order 100 than that needed for thermal 
freeze-out. 

A distinctive feature of the model is that small mass 
splittings SM ~ 1 MeV between the DM states are gen- 
erated at one loop. Such splittings can be important 
for making the model consistent with direct detection 
constraints, by causing DM scattering on nucleons to be 
inelastic and endothermic. This allows the rate of DM 
interactions in detectors to be sufficiently suppressed de- 
spite having relatively large couplings to nucleons 33]. 
The mass splittings can moreover lead to a significant in- 
crease in the Sommerfeld enhancement factor 38] , which 
effect we take into account here. In addition, they can po- 
tentially help to explain the annual modulation observed 
by DAMA/LIBRA through the inelastic DM mechanism 
(see for example [39J ) , or the 511 keV excess observed by 
INTEGRAL j|o| using the excited DM mechanism [Zlj| - 
44], although we do not explore these directions in the 
present work. 

A further refinement we make compared to most pre- 
vious analyses is to take into account the dependence 
of the Sommerfeld enhancement on the position in the 



Galaxy, an effect that was first pointed out in ref. [45j. 
The Sommerfeld enhancement depends sensitively on the 
DM velocity, and both the DM velocity dispersion and 
escape velocity depend upon r. We find that this typ- 
ically makes the enhancement larger near the galactic 
center where the gamma rays that are most important 
for the IC constraint are produced. The constraint thus 
becomes harder to satisfy than in models with a spatially 
constant boost factor, as is usually assumed. 

Nevertheless, we find examples of models that are 
marginally able to produce leptons consistent with the 
PAMELA and Fermi observations, and which are also 
barely consistent with the various gamma ray constraints. 
The preferred models have DM with M = 800 GeV and 
gauge boson masses below 200 MeV, so that only elec- 
trons and no muons or charged pions are produced in 
the annihilations. However we find that the Sommerfeld 
enhancement factor would always be too large to achieve 
this concordance if the DM was present at the expected 
relic density. We are obliged to suppress the rate of anni- 
hilations by introducing a factor 1/ / in the density of the 
annihilating component. This can be achieved by increas- 
ing the gauge coupling appropriately so that the relic 
density is reduced. The minimum reduction 1/ / = 0.4 
occurs if the DM mass splitting is negligible, and smaller 
values 1// =* 0.14 - 0.2 are needed if SM ~ 1 MeV. 
This comes about because of the JM-dependence of the 
Sommerfeld enhancement factor for multistate DM anni- 
hilations [381 ] - 

In the present work we do not consider astrophysical 
boost factors due to increased annihilation rate in sub- 
structures of the main halo [46|-[52j. Inclusion of these 
effects would presumably strengthen the gamma ray con- 
straints on the models, but on the other hand a fully 
consistent treatment of substructure effects should simul- 
taneously consider the modification to the lepton signal 
from annihilations in subhalos [5J|. This could conceiv- 
ably have a compensating effect that would weaken the 
constraints due to requiring less e production near the 
galactic center [54{ . Further investigation of this issue is 
in progress. 

We start by introducing the particle physics models 
in section [Til and recalling the constraint upon the dark 
SU(2) gauge coupling imposed by the relic density. In 
section IIIII we analyze how large an enhancement factor 
is needed for the models to explain the anomalous lepton 
observations, the theoretical computation of the Som- 
merfeld enhancement, and its phase space average in the 
galaxy. Section IIVI explains how we implement the in- 
verse Compton 7 ray constraint for general DM Einasto 
profiles, while other relevant constraints are discussed in 
section [V] Our main results showing which parameter 
ranges can be compatible with all the observations are 
presented in section |VT1 We discuss their implications in 
section lVTll Three appendices give details for the branch- 
ing ratio of ground state DM annihilation into leptons 
(app.[X]), the multistate Sommerfeld enhancement (app. 
IB1 and the calculation of the DM escape velocity (app. 



[C|l as a function of radial position. 



II. PARTICLE PHYSICS MODELS 

A. The models 

For the purpose of testing a concrete and specific 
model, we will analyze the simplest examples of non- 
abelian hidden sector models that can explain the var- 
ious cosmic ray anomalies; namely a dark SU(2) gauge 
group under which the DM transforms as a doublet or 
a triplet [35l [36| . The Lagrangian for the DM and the 
gauge bosons is 

c = \H0ij - MMxj - ^b; v b^ - i A„sr v 

(i) 

We assume that the dark Higgs triplet A a gets a VEV in 
the a = 2 direction, so that B2 mixes with the standard 
model hypercharge Y^ v with strength 



= (A 2 )/A 



(2) 



This causes B2 to acquire a coupling of strength eq to any 
SM particle of charge q, and therefore to mediate the de- 
cay %2 — ► Xi e+e ~ m the doublet model, or X3 — ► Xi e+e ~ 
in the triplet model. 1 The value of e is not strongly con- 
strained; it can be anywhere between 10 -2 and 10 -9 [36[ 
(except for a small region around 10 -6 -10 -7 excluded by 
the E137 experiment if the dark gauge boson mass is less 
than 400 MeV [H). 

These models can explain the PAMELA positron ex- 
cess and Fermi/LAT e^ excess through the annihilation 
X1X1 — > B2B2, followed by the decays B2 —> ff, where 
/ is any charged SM particle with mass less than /12/2 
(where P2 denotes the mass of B^)- Although higher 
mass fermions can be produced through off-shell gauge 
bosons, such processes have a much smaller cross section 
than those in which the B^s are on shell. Using this fact, 
one can account for the lack of any antiproton excess in 
the PAMELA observations by assuming that fj,?. is less 
than 2m„. 



corresponding to g = 0.98, 0.62 respectively at M = 1 
TeV. We will adopt this relation between g and M in our 
analysis of the triplet model. For the doublet model, the 
DM must be Dirac in order to have a bare M\X mass 
term that respects the SU(2) gauge symmetry. It could 
therefore have a chemical potential that fixes its relic den- 
sity independently of its annihilation cross section. This 
could provide motivation to consider a g as an additional 
free parameter for the doublet model, if one can find a 
mechanism for producing the asymmetry between x and 

X- 

In later sections we will find that it is difficult to satisfy 
all constraints if the annihilating DM has the full allowed 
relic density, because the Sommerfeld enhancement of the 
annihilation cross section is too large. A more general 
possibility is that the gauge coupling exceeds ^ by a 
factor y/J, which leads to an enhancement in the cross 
section a — > fa, and a consequent suppression in the 
density by p — > p/f. Since the signals scale like p 2 a, this 
leads to a reduction by I// in the rate of annihilations in 
the galaxy. Thus more generally we will consider gauge 
couplings given by 



V7 



M 



1 TeV 



0.077, doublet 
0.031, triplet 



(4) 



where / > 1. 

If / > 2, then less than half of the DM resides in 
the ground state, and there must exist some additional 
DM sector to make up the rest of the mass density. If 
I < / < 2, the needed reduction could happen naturally 
within the model due to the fact that some fraction (less 
than 1/2) of the DM resides in one of the excited states, 
assuming it is stable, and that it is not able to annihilate 
into leptons. (Our triplet DM model provides such an 
example.) A quantitative calculation of the relic density 
of the excited state is beyond the scope of this paper, but 
will be carried out in ref. |44|. 



C. Higgs sector 



B. Relic density 

Ref. [36| computed the cross section for XX ~^ BB and 
determined that it leads to the right relic density from 
thermal freezeout in the doublet or triplet model, if the 
gauge coupling a g = g 2 /An has strength 



—) 

1 TeV/ 



0.077, doublet 
0.031, triplet 



(3) 



The gauge couplings have the form 9Xi a ti^oXi or 
29 e abcXa$bXc respectively, in the doublet and triplet models. 



The model must also include a dark Higgs sector to 
break the SU(2) gauge symmetry and give masses p £ 1 
GeV to the gauge bosons. One-loop self-energy correc- 
tions to the DM masses induce mass splittings of or- 
der SM ~ a g 8p if the gauge boson masses themselves 
are split by an amount of order 8 p. The splittings 
SM ~ 1 MeV are potentially useful for explaining the 
INTEGRAL 511 keV gamma rays by the excited dark 
matter mechanism, and they are needed to satisfy con- 
straints from direct detection unless the mixing parame- 
ter e defined in ^ is small, e < 10~ 6 [33]. Larger mixing 
e ~ 10~ 3 are more interesting from the perspective of 
laboratory tests by fixed-target experiments. 

In the simplest and most predictive models, the DM 
states would have no Yukawa couplings to Higgs bosons 



of the hidden sector. 2 If one is also interested in trying 
to address the 511 keV gamma ray excess, the relation 
SM = ^a g Sfi gives a constraint on splittings between 
the gauge boson masses. The DM mass splittings should 
not be much more than of order 1 MeV; otherise excited 
DM decays would produce e of too high energy to be 
compatible with the narrow observed 511 keV spectrum. 
For the triplet DM model with a g = 0.031 for M = 1 
TeV, we thus expect that <5/x £ 180 MeV if the largest 
DM mass splitting is SM — 3 MeV. Moreover we will 
show that models with /i ~ 200 MeV can more easily 
satisfy stringent 7 ray constraints. 



D. Explicit example 

Ref. [36| showed that the minimal Higgs sector for 
splitting triplet DM is one that gives a VEV to an ad- 
ditional triplet Higgs field A' a in a direction orthogonal 
to that of the first triplet VEV, e.g., A[. We do not 
wish to commit ourselves to a particular Higgs potential 
in the present analysis. Instead, we will assume that it 
is possible to design one that gives the desired DM mass 
splittings of order 100 keV to a few MeV, and gauge bo- 
son masses of order a few 100 MeV. In fact, this is quite 
reasonable even in the context of the simplest model with 
two Higgs triplets. Using the analysis of ref. [3g|, the 
mass splittings SMi 3 = 3.1 MeV, SM 23 = 100 keV and 
gauge coupling g = 0.62 imply that Aj = 320 MeV, 
A2 = 80 MeV, and the gauge boson masses are /Ui = 50 
MeV, [12 — H3 — 200 MeV using this simplest Higgs sec- 
tor. 



III. SOMMERFELD ENHANCEMENT FOR 
PAMELA/FERMI LEPTON SIGNALS 

In this section we explain how we calculate the effective 
cross section for annihilating DM to produce high-energy 
e (and possibly heavier charged particles) in the galaxy, 
relevant for explaining the PAMELA and Fermi e ± ex- 
cesses. A particularly simple case to focus upon is that 
where the mass /i of the B^ gauge boson is below the 
threshold for producing a [i + -\x~ pair. In that case one 
need only consider final states consisting of e . Moreover 
the cross section to fit the PAMELA/Fermi lepton data 
is smaller in this case, and thus easier to achieve from 
the model-building perspective. Another reason for pre- 
ferring the 4e channel is that it is not ruled out by CMB 
reionization constraints, whereas the 4/i channel is (de- 
pending upon the assumed value of the local DM density 
Po, as we will discuss); this is apparent in figure [T] 
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FIG. 1: Heavy lines: allowed regions for explaining the 
PAMELA and Fermi lepton excesses, assuming 4e or i/j, final 
states, from ref. [171 ] . The fits assume an Einasto DM profile 
with p© = 0.3 GeV/cm 3 , a = 0.17 and r s = 20 kpc. Narrow 



lines show constraints from CMB [24 
gamma rays [171 ] . 



and inverse Compton 



A. Cross section ratio for X1X1 — ► e 

Although the models we consider could have an ex- 
cited DM state that is stable or metastable, the ground 
state normally should have a larger density m due to 
downscattering processes like X2X2 ~^ X1X1 m the early 
universe. To simplify the analysis, we assume that it is 
a good approximation to ignore annihilations involving 
the subdominant component for the production of lep- 
tons via \\ — > 4e. For the triplet DM model, there are 
in fact no such annihilations in the form X2X2 —> 4e be- 
cause xi does not couple to B2, the gauge boson that we 
have assumed to be the only one mixing with the SM. 
The only annihilations we miss in this case are of the 
form X1X2 - ► B\e + e~ , where the B\ appears as missing 
energy. 

It is potentially important to realize that only the 
ground state (and possibly some small addition of ex- 
cited state particles) contribute to the annihilation in the 
galaxy, because the cross section for X1X1 to annihilate is 
different from that for XiXj m the early universe. They 
are related by the ratio &n, 



0-(XiXi -► e ) = buvfaiXj -> an y) 



which in the doublet and triplet models is given by 

|, doublet 

|, triplet 
as we explain in appendix [Al 3 



611 = 



(5) 



(6) 



2 In a model with only triplet Higgses and Majorana DM, the only 
potentially gauge invariant coupling £ a bcXa.Xb^c in f ac t vanishes. 
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FIG. 2: Upper left: Sommerfeld enhancement factor (|B1[) as a function of gauge boson mass fi for DM mass M = 1 TeV, 
mass splitting 8M — 1 MeV, and gauge coupling a g — 0.031. Different curves are labeled by DM velocity in center of mass 
frame, in km/s. Upper right: the same, but averaged over the DM velocity distribution functions as in (jlOp . for three different 
circular velocities vq — 200, 240 and 280 km/s. The allowed regions for Fermi with 4e annihilation channel for the doublet and 
triplet DM models, assuming solar density pq = 0.3 GeV/cm and DM fraction / = 1 (see eq. Q) are shown by horizontal 
lines. Central panels: same as upper but with vanishing mass splitting 8M = 0. Lower panels: velocity-averaged Sommerfeld 
enhancement factors versus fi for doublet model, with a g = 0.077 and SM = 1 MeV (left), 8M — (right). 



computation of this quantity 



B. Sommerfeld enhancement for multistate DM 

The boost factor needed for explaining the 
PAMELA/Fermi e^ observations can be inferred 
using figure [TJ which shows the 3er-allowed regions 
for M and av, assuming 4e and Afi final states, and 



also assuming the solar neighborhood DM density to 
be Pq = 0.3 GeV/cm 3 . The central value of the 4e 
allowed region has av — io~ 231±0 - 2 cm 3 /s at M = 1.1 
TeV, while the av needed for the right relic abundance 
is 3 x 1CP 26 cm 3 /s. This leads to a required boost 
factor of B = l(T 231±0 - 2 /(>ii • 3 x 1CT 26 ). Allowing for 
av — > fav and p — > p/f as described in section III Bl the 
required boost factor scales as B — > fB. For the doublet 
and triplet models, using the gauge couplings ((4]), this 
gives 



log B 4e = log / 



2.67 ±0.2, doublet 
2.36 ±0.2, triplet 



(7) 



for the boost factor needed in the 4e annihilation channel. 

To get a feeling for the ability of the models to achieve 
large enough boosts, we start by computing the Sommer- 
feld enhancement factor for a fixed DM mass of M = 1 
TeV near the central region allowed by Fermi with 4e 
final states. For this purpose one should take into ac- 
count that there is more than one DM state, so the usual 
Sommerfeld factor for single-state DM is not appropriate. 
The general case has not been solved explicitly, but that 
of two states with therefore a single mass splitting was 
recently analyzed in [38j , where it was found that the en- 
hancement factor can be approximated by eq. (|B1[) given 
in appendix [Bl 

The result (jBll) applies directly to the doublet model, 
where there are only two states, but only approximately 
in the case of the triplet model, since the latter has three 
states and two mass splittings. Since the enhancement 
is an increasing function of SM, one might reasonably 
expect the multistate effect in the triplet model to be 
captured by choosing the largest of the two mass split- 
tings. We make this assumption, which in a more ambi- 
tious study should be checked. It is however a technically 
challenging problem. 



C. Velocity-averaged enhancement at r — r© 



in the region r < 10 kpc. (We elaborate on the r depen- 
dence of vo in section fVI At ) At the solar radius, which 
we take to be r Q = 8.33 kpc [58|], this gives v csc — 2.53 vo, 
corresponding to escape velocities in the range 440 — 610 
km/s. This is in reasonable agreement with constraints 
from direct measurements of high- velocity stars [59| . We 
need to average S(\vi — #2 1) over the phase space, 

(S) = /rfVd 3 v 2 /(fi)/(« 2 )^(|t/i -tfc|) (10) 



D. Predicted versus desired enhancement factor 
values 



In figurc[2]the dependence of the Sommerfeld enhance- 
ment S and its velocity average (S) on the gauge boson 
mass p is shown, for the example of M — 1 TeV and for 
the triplet model with a g = 0.031, as well as the doublet 
with a g = 0.077. These are the gauge couplings needed 
for the right relic density in the respective cases. The ex- 
tra enhancement due to the mass splitting SM — 1 MeV 
is quite significant; in fact it tends to give rise to boost 
factors that are too large compared to the values needed 
to explain the high-energy lepton excesses. The effect 
of SM actually saturates around SM = 600 keV for the 
triplet model and SM = 100 keV for the doublet model; 
the enhancement remains roughly constant for larger val- 
ues (however the validity of the approximations leading 
to JET) breaks down if SM > a 2 g M). 

We conclude that, from the point of view getting the 
right boost factor, there is no need for the mass splitting, 
and in fact too large a mass splitting is disfavored. How- 
ever this can be compensated by invoking the /-factor of 
section ITTB1 since taking / > 1 increases the value of the 
needed boost factor. 



We must compute the boost factor in the neighbor- 
hood of the Sun in order to predict the rate of high- 
energy lepton production from DM annihilations. The 
Sommerfeld enhancement depends upon the relative ve- 
locity of the DM particles, which is usually characterized 
by a Maxwellian distribution 

f( v ) = Ne- v2 / v ° 9{v - v csc ) (8) 

with a cutoff for v above some escape velocity v csc . The 
value vo is commonly taken to be 220 km/s at the so- 
lar radius r Q = 8.3 kpc, although higher values ~ 250 
km/s have been advocated more recently [56| . The range 
wo € 200 — 280 km/s is suggested by the compilation of 
different measurements in ref. 57|. The escape velocity 
is directly correlated with vq , as detailed in appendix [D] 
There we motivate our choice for the r-dependent rela- 
tion 



z; 2 sc (r) = 2vl{r) [2.39 + ln(10 kpc/r)] 



(9) 



IV. INVERSE COMPTON GAMMA RAY 
CONSTRAINT 



To undertake a thorough exploration of the parameter 
space, we must take into account constraints that could 
rule out the models. The intepretation of the PAMELA 
and Fermi excesses as being due to DM annihilation has 
come under pressure from a number of complementary 
constraints. Many of these arise from radiation that 
would be produced by the leptons after annihilation, 
which should be directly detectable as gamma rays or 
radio emission in the galaxy, or indirectly in distortions 
of the cosmic microwave background by contributing to 
the reionization of the early universe. The most stringent 
constraint on our model is from inverse Compton scat- 
tering in the galaxy. We therefore consider it separately 
in this section, and discuss other constraints in section 

El 



IC bound on en>(DM DM -» /i + /T) in 10~ 23 cm 3 /sec for M = 1.3 TeV 
EINASTO DM profile with L = 4 kpc 
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FIG. 3: Regions near galactic center used by ref. [171 ] for computing the Fermi IC bound. Numbers in each region refer to 
constraints arising for a particular annihilation channel (into /i + /i_), but size and shape of regions are the same for the channel 
we consider, 2e + 2e~. Courtesy of A. Strumia. 



Characterization of the bound 



If DM annihilates into charged particles, the latter will 
undergo inverse Compton scattering on the galactic radi- 
ation field (for example, starlight) to produce high energy 
gamma rays that could be detected by the Fermi LAT. 
Several authors have considered the constraints arising 
from the nondetection of such a signal (13T - I20J . Most of 
these assume a two-lepton final state rather than four 
leptons, as is the case in our models. The constraints on 
the former are stronger because their spectrum is harder 
than in the four-body case. Ref. [TtJ has considered four- 
lepton final states for several DM density profiles, so we 
adopt their results for the present analysis. 

The resulting fits of P3 for the PAMELA/Fermi ex- 
cesses and the IC constraints for the 4e and 4p channels 
are shown in figure [TJ for a particular choice of Einasto 
profile, 



p(r) 



p exp 



(11) 



6 
5 




i 

~* *"• ^ 

1 




- 


n 4 


r s = 10 kpc 

__ / = 30kpc 


- 


Q_ 
Q. 3 

O 

- 1 2 

1 





^v 



-0.5 

log r (kpc) 



0.5 



FIG. 4: Einasto density profiles normalized to pq for a — 
0.17, 0.20, 0.23 and r s = 10 kpc (solid curves) or r s — 30 kpc 
(dashed curves) 



section for these two cases can be approximated as 4 

4e:loger 4e v = -15.14 - 14.66 x + 8.00 x 2 - 1.78 x 

+ 0.144 a; 4 
4/i : log a^v = 



245 - 21.37 a; + 10.43 x 2 - 2.136 a; 3 



0.161a; 4 



(12) 



where av is in units of cm 3 /s and x — logM with M 



with p G = 0.3 GeV/cm 3 , a = 0.17, r s = 20 kpc, and 
r Q = 8.5 kpc. The maximum allowed annihilation cross 



tables on which these fits are based provided by A. Strumia 
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FIG. 5: Contours of M6o/(10 11 Mq) (Meo being the mass contained within 60 kpc) in the plane of a and r s (in kpc units), for 
several values of p©. The heavy line labeled "4" is the central value, and other heavy lines (at 4 ± 0.7 and 4 ± 1.4) denote the 
neighboring la and 2<r confidence intervals determined by ref. [701 ]. 



in GeV. The constraints are found by comparing the ob- 
served gamma ray spectrum from Fermi data with the 
predicted IC contribution from the DM-produced lep- 
tons, and demanding that the latter do not exceed the 



former in any energy bins by more than 3a. Ref. 17 1 
carried out this procedure for different regions of the sky, 
similar to those depicted in fig. |3] The regions closest to 
the galactic center are the ones that give the strongest 
constraint on the cross section, resulting in the bounds 
(fT2|) . It is worthwhile to notice that the inner 5° of 
galactic latitude have not been used in deriving the con- 
straints, due to the difficulties of interpreting foregrounds 
associated with the disk of the galaxy. We have adopted 
ref. [13] 's constraints relative to the 'MED' propagation 
model for e ± , corresponding to a diffusion zone thickness 
of L = 4 kpc; the results from other choices (except for 
the extreme value L — 1 kpc) are not significantly differ- 
ent in the range of DM masses in which we arc interested. 

For realistic models of the type we consider, if the 
gauge boson is heavy enough to decay into muons, then 
electrons will also necessarily be produced with some 
nonnegligible probability, depending only on the relative 
phase spaces for decays into the two different final states. 
Moreover charged pions will also be produced if the gauge 
boson mass /u is greater than 2to w . In such cases the ac- 
tual bound will be somewhere in between cr 4e and a^ 
given in (fT2")l . We will avoid the uncertainty of how to 
correctly interpolate between the two bounds by restrict- 



ing our attention to the case /i < 2m M , so that only the 
4e bound applies. 



B. Dependence of IC bound on DM profile 

Figure Q] shows that for the assumed fiducial density 
profile (TTTJ) , the constraints rule out the Fermi-allowed 
regions. This can be circumvented however by consider- 
ing DM density profiles that are less strongly peaked near 
the galactic center, where the most stringent constraints 
come from. With the Einasto profile, both a and r s have 
an effect on the cuspiness of p(r) near the galactic center, 
which is illustrated in figure [4j larger values of a or r s 
weaken the bound. One of our goals is to quantify this 
statement to determine the range of a and r s that gives a 
consistent description of the PAMELA/Fermi anomalies 
in terms of DM annihilation. 



1. Theoretical and observational constraints on the DM 
profile 

Let us first consider what are the reasonable ranges of 
variation for the halo profile parameters a and r s based 
on theoretical and observational considerations. iV-body 
simulations like Aquarius find 0.115 < a < 0.179 [32|]; in 
a different fit of the same galaxies over a slightly smaller 
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FIG. 6: Representation of the quantites r, s, ip used to com- 
pute the J factor (|14[) . The coordinate s runs along the line 
of sight in the direction specified by the angle ip and is related 
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radial region, ref. [60| finds a as large as 0.19. These 
values are for pure DM simulations without baryons; the 
inclusion of baryons tends to reduce a dramatically to the 
range 0.07 — 0.015 due to the concentrating effect of the 
baryons on the inner halo |60h63 |. In contrast, the pure 
DM simulations find a larger range of r s , from 15 kpc to 
29 kpc (where we have assumed h = 0.7 for the Hubble 
parameter, since the results are quoted in units of kpc/ ft.). 
This range also gets decreased when baryons are in- 
cluded, to 5— 15 kpc. There thus appears to be a conflict 
between the DM simulations with baryons (BDM) and 
the annihilating DM interpretation of PAMELA/Fermi 
events, since the adoption of a — 0.15 and r s = 15 kpc 
will only make the IC constraint stronger, compared to 
values that already rule out the model. We also note that 
one of the highest resolution simulations, GHALO [65| . 
obtains a = 0.155 even without including baryons. 

For the allowed range of the solar neighborhood den- 
sity pq , observations provide tighter constraints than do 
simulations. The value p©^ = 0.3 GeV/cm 3 has for a 
long time been considered standard [66J]. More recently 
a higher central value 0.385 ± 0.027 GeV/cm 3 has been 
advocated in [67( and 0.43 ± 0.15 GeV/cm 3 has been de- 
termined using a method that does not rely upon the 
detailed form of p(r) in |68J . 

Observational limits on the ranges of the other Einasto 
parameters tend to be weaker than the ranges suggested 
by simulations, but this depends upon the assumed value 
of pq . One constraint comes from the mass within a given 
radius inferred from rotation curves of the Milky Way 
galaxy. Ref. [69| inferred a total mass within 50 kpc of 
the galactic center, M50 = 5.3 x 10 11 M©. More recently, 
using data from the Sloan Digital Sky Survey, ref. pZQj ] 
obtained the mass within 60 kpc as Mqq — (4 ± 0.7) x 
10 11 M Q . In fig. [5] we plot contours of Af 60 /(10 n Af Q ) 
in the a-r s plane, for several values of pq. At p© = 0.4 
GeV/cm 3 , values of r s exceeding 20 kpc are disfavored by 
this measurement, while at p© = 0.3 GeV/cm 3 it imposes 
essentially no constraint. 



FIG. 7: J factor eq. ()14[) over the fiducial value used by ref. 
[17J as a function of a, for several values of r s . Dashed curves 
are the analytic approximation (|15[1 . 



2. J factors for different density profiles 

A principal goal of this work is to investigate the sen- 
sitivity of the IC bounds in (fT2")l to different choices of 
the Einasto profile parameters. The dependence upon 
Pq is easy to quantify, since the rate of DM annihilations 
is proportional to p 2 a. This implies that for a general 
value of pq, and the fraction 1// of annihilating DM de- 
fined through (j4|, the upper bound on the annihilation 
cross section scales as 



0"4e 



C4e 



PO,0 
P& 



f 



(13) 



relative to that at the reference solar density p 0! o = 0.3 
GeV/cm 3 and / = 1. This means that the constraints 
become more severe by a factor of (0.43/0.3) 2 = 2 if 
we adopt the new central value 0.43 GeV/cm 3 of ref. 

68 1 instead of the value 0.3 GeV/cm 3 assumed in ref. 

17]. However one must remember that the cross section 
required to fit the PAMELA and Fermi leptons goes down 
by the same factor, so in fact the value of p© has little 
direct effect on the ability of a model to satisfy the IC 
constraint. 5 

The effect of changing the shape of the profile, on the 
other hand, is harder to quantify. The correct procedure 
would require solving the diffusion equation for electrons 
in the galaxy again for each choice of a and r s , and then 
recomputing the spectrum of IC scattering. Here we con- 
tent ourselves with a simpler method of estimation, in 
which we ignore the changes to the electron diffusion in- 
duced by changes in the shape of the profile. In this 
approximation, the rate of IC scattering due to leptons 



Indirectly it has an effect through the constraint from Mqo , which 
makes large values of r s hard to achieve when pQ is large. 
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originating from the DM annihilations changes only by a 
geometrical factor, the averaged J factor, defined as 



J 



1 



ds 

l.o.s. r O 



dn / — 



p© 



(14) 



where the integral along the line of sight is then averaged 
in rp (defined in fig. |6]) over the solid angular regions Ail 
that are used to compute the gamma ray constraints, as 
shown in fig. |3l As noted above, the relevant regions 
for obtaining the strongest bound are those with latitude 
5° < |6| < 10° and longitude \l\ < 10°. The direction 
tp corresponds to latitude b and longitude £ as cos?/' = 
cos b cos i. 

The resulting J factors, divided by the fiducial one Jo 
corresponding to the parameters chosen by ref. [13|, are 
shown as a function of a for several values of r s in fig. 
[7J The IC upper bound on the annihilation cross sec- 
tion is expected to scale like o± e — > cr ie /{J / 'Jo), so that 
smaller values of J/ Jo correspond to a weaker bound. 
One observes that it is not easy to weaken the bound 
while remaining within the expectations of iV-body sim- 
ulations. To achieve a factor of 2 reduction would require 
a ~ 0.22 and r s = 30, for example, which appear to be 
extreme values. 

One can get analytic insight into the dependence of J 
on a and r s by hypothesizing that J scales like 



_£_ ^ / po(nc) V 
Jo V p(nc) / 



(15) 



where po is the Einasto profile using the fiducial a = 0.17 
and r s = 20 kpc values, and rjc is some characteristic 
radius that should be of order r Q times the angular dis- 
placement (i.e., 5°) from the galactic center of the rel- 
evant solid angular region. By tuning nc to the value 
1.75 kpc, we are able to get good agreement with the 
numerical results, with less than 4% error in the region 
a > 0.15 (and less than 6% elsewhere). The approxima- 
tion is shown as the dashed curves in figure [7) 



3. Potential for strenghtemng of IC constraint 

As we noted above, the current constraints are based 
on a region that excludes the central latitudes |6| < 5°, 
due to the complexity of the disk of the galaxy. On the 
other hand we know that the strongest bounds arise from 
the innermost regions. It is therefore interesting to try to 
project how much stronger the constraints might become 
if one had made use of the Fermi/LAT 7 ray data from 
the inner latitudes. It is straightforward to recompute 
the J factors over regions that include the inner ±5° of 
latitude to try to project how much stronger the bounds 
might become. We show the result in fig. [3 According 
to this extrapolation, the bound could get stronger by a 
factor of 2 at the reference values of a = 0.17 and r s = 20 
kpc. However this is merely suggestive, and one would 
have to analyze the actual Fermi 7-ray data in the inner 
region to draw firm conclusions. 




FIG. 8: Solid curves: J factors as in fig. [7] but averaged over 
the angular region including the inner latitudes of the galaxy. 
Dashed curves are for the original region considered in figure 
[7J for comparison. 



V. OTHER CONSTRAINTS 

As mentioned above, a number of complementary con- 
straints have been studied pertaining to models of annihi- 
lating DM. We present here the most significant of them, 
and discuss how they apply to the cases in which we are 
interested and their impact on the parameter space. 



A. GC gamma ray constraints 

Fluxes of gamma rays are directly produced by the 
DM annihilation process itself ('prompt' gamma rays), 
mainly from the bremsstrahlung of charged particles and 
the fragmentation of hadrons, e.g., tt°, produced in the 
annihilations. These fluxes extend to very high energy 
(up to the mass of the annihilating DM particle) and 
are therefore better constrained by comparing the pre- 
dicted flux with the observations of high energy gamma 
ray telescopes such as the H.E.S.S. or VERITAS (Very 
Energetic Radiation Imaging Telescope Array System) 
observatories. Refs. @ and [8|, [lTj have indeed carried 
out the analysis of these constraints on the parameter 
space of models that aim to explain the lepton anoma- 
lies. H.E.S.S. observations of the region of the Galactic 
Center [71| (defined as a disk of radius 0.1° centered at 
the GC) , of the Galactic Ridge 72] (defined as a region 
of 0.8° x 0.3° in longitude and latitude centered at the 
GC) and of satellite dwarf spheroidal galaxies such as 
Sagittarius, Draco, Ursa Minor or Willman 1 17311 impose 
constraints that can be significant (see [8|, Mi LL2I ) ■ 

However, all of these constraints are very sensitive to 
the details of the chosen DM profile. In particular, for 
the Milky Way ones, it is the inner part of the galactic 
DM halo which counts, and numerical simulations do not 
provide direct determinations of the profile at r ~ 100 
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pc. Moreover, assuming that an extrapolation is possi- 
ble, only for the steepest of the profiles that we consider 
(small a values) do these bounds become competitive 
with those from IC, and such values are anyway disfa- 
vored by the more robust IC constraints. We therefore 
need not consider these bounds any further. 



B. GC radio constraint 

The e produced by DM annihilations within the 
galactic magnetic field emit synchrotron radiation, which 
falls in a range of frequencies roughly spanning the ra- 
dio to the IR. The Galactic Center is the best region to 
search for this effect, both because of the large local value 
of the DM density and of the magnetic fields. We do not 
enter here into the details of the needed astrophysical and 
particle physics ingredients but refer to [9( for the com- 
plete discussion. There, in brief, the signal is computed 
neglecting advection and diffusion but scanning different 
(extremal) assumptions for the galactic magnetic field. 
Ref. [8] has in particular considered the case of the 4- 
lepton annihilation final state that we interested in. 

Comparing the predicted flux with observations pro- 
duces constraints on the DM annihilation cross section. 
Since the observed GC microwave spectrum is harder 
than what DM annihilations can produce, the dominant 
bound is obtained considering the observation available 
at the lowest observed frequency, v — 0.408 GHz, per- 
formed by [74j | in a region with full width half maxi- 
mum of 4", corresponding to about 0.1 pc. The resulting 
bounds can be quite stringent (see [8[ ) and the constraint 
extends to low DM masses where the 7-ray bounds from 
H.E.S.S., discussed above, are not effective. The varia- 
tion of the magnetic field negligibly affects the bound, 
because the radio emission is predominantly produced 
by outer regions. A subdominant bound comes from 
the VLT observation [75| at the larger infrared/visible 
frequency from a region with even smaller angular size 
0.04" i.e., r < 0.0016 pc. This bound somewhat depends 
on the magnetic field profile, and it becomes numerically 
significant only for spiked DM density profiles [76j ]. 

However the same discussion as above applies to these 
bounds: as for the GC gamma ray constraints, these 
bounds are very sensitive to the details of the chosen DM 
profile (they originate from the even smaller regions of 
r ~ 1 pc around the GC) and do not apply to non-steep 
profiles. We therefore need not consider these bounds 
any further. 



photons directly emitted in the annihilation of two DM 
particles, and by the lower energy photons produced by 
inverse Compton scattering. The latter turns out to be 
by far the most important process; in fact, the cross sec- 
tion for je~ scattering decreases rapidly with the en- 
ergy of the impinging photon, so that low energy pho- 
tons are more efficient in removing the electrons from 
the atoms. These 'primary reionization' electrons then 
deposit their energy in the intergalactic medium via sev- 
eral other interactions, freeing many more electrons and 
also augmenting the temperature of the gas. One way 
to constrain DM annihilation properties is therefore to 
look at the modifications of the CMB spectrum pro- 
duced by the fact that the CMB photons meet an opaque 
medium (with free electrons produced by the ionization) 
in their journey from the surface of last scattering. This 
is in particular encoded in the total optical depth of the 
Universe r parameter, r is measured by WMAP to be 
t = 0.084 ± 0.016 [77], of which about 0.038 is due to 
the low-redshift reionization (z < 6) produced by stars. 
A DM-induced optical depth larger than 0.062 (the la 
upper bound of [77| ) is therefore excluded by these argu- 
ments. 

This bound has been studied in the literature in several 
references [2l| - [25| . with mutually consistent results. The 
fit to the numerical result in [24| is 



logfTcMBf = —26.3 + 1.15 log M 



(16) 



where av is in cm 3 /s and M is in GeV. While this was 
originally computed for 2-lepton final states, it applies 
equally well for the 4-lepton ones in which we are in- 
terested. All that matters is the total amount of en- 
ergy that is injected in the primordial intergalactic gas 
in the form of "electromagnetically coupled" final prod- 
ucts (electrons, positrons, photons, possibly hadrons, but 
not neutrinos). All such final products end up causing 
reionization and thus producing free electrons, regard- 
less of the precise way in which the annihilation occurs. 
We neglect minor differences possibly introduced by the 
difference in the shape of the spectra. The dotted lines 
in fig. Q] and QT] show these bounds. What one sees is 
that, for the 4e case in which we are interested, they rule 
out the large DM mass portion of the PAMELA-only fit 
region, but not the PAMELA+Fermi one. 

In the present work, we consider the case where the 
annihilating DM constitutes a fraction 1// of the total 
DM. Since the annihilation rate scales like p 2 <r, the bound 
(fTBf will be weakened by ctcmb — > / 2 ccmb- 



C. Reionization constraint from CMB 

The flux of energy injected by DM annihilation, from 
the recombination epoch until today through the forma- 
tion history of DM halos, results in ionization and heating 
of the intergalactic medium. The ionization and heating 
can be produced both by the highly energetic 'prompt' 



D. Extragalactic gamma ray constraints 

Finally, we briefly review the extragalactic gamma ray 
constraints on DM annilation. These refer to the fluxes 
of prompt and inverse Compton gamma rays that are 
produced by annihilations in all DM halos outside of our 
own galaxy and throughout the history of cosmological 
structure formation. They reach us (properly redshifted) 
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FIG. 9: (a) Left: circular velocity (|17(l as a function of r for exponents \ = 1-64 (includes effects of baryons) and % = 1.95 (pure 
dark matter) at standard Einasto parameters a — 0.17 and r a — 20 kpc, and for a very cuspy profile (dot-dashed curve), (b) 
Right : velocity-averaged Sommerfeld enhancement as a function of r for two different values of the circular velocity v c = vo (r© ) . 
The ellipses denote regions in r where the dominantly observed IC 7 rays and respectively leptons are produced. 



uniformly from all directions, so that they have to be 
compared with the Fermi measurements of isotropic dif- 
fuse gamma rays [78| . A number of references have gone 
through this analysis [lj|, [20, UQ, [27[ , showing that the 
bounds depend very significantly on assumptions about 
the parameters of the history of structure formation. 
E.g., for conservative choices of the halo concentration 
parameter function, all of the parameter space in which 
we are interested is allowed by these bounds. 



VI. RECONCILING ANOMALOUS LEPTONS 
WITH 7-RAY CONSTRAINTS 

We have now discussed most of the necessary ingredi- 
ents for scanning over the parameters of the theoretical 
models and of the DM density and velocity distributions, 
to see if there exist any choices that are consistent with 
all the observational requirements. It remains to explain 
some details of our methodology for predicting the ef- 
fective annihilation cross sections, both for the lepton 
signal and for the IC 7 ray signal. These can differ from 
each other somewhat because of the fact that most of the 
observed leptons are produced within 1 kpc of the Sun, 
whereas the most important IC 7 rays come from closer 
to the galactic center, and the Sommerfeld enhancement 
of the cross section is position dependent. We explain 
these details in the present section, and use them to ob- 
tain the final results. 



A. Position dependence of enhancement factor 



As explained above eq. (jlO)) , the Sommerfeld enhance- 
ment must be averaged over the phase space of the DM. 
Previously we only considered this in the vicinity of the 
Sun, but now we must take into account that the average 



DM velocity as well as the escape velocity depend upon r. 
This dependence has been measured in iV-body simula- 
tions, and is shown to follow a scaling relation predicted 
by Bertschinger [791 ] . 



vo(r) 3 



' p(r) 



(17) 



where the exponent \ = 1.9 — 2.0 for pure DM simula- 
tions [32], and it takes lower values \ = 1.60 — 1.67 in 
simulations including baryons [60J. The r-dependence of 
the escape velocity is estimated as in eq. ©■ We plot the 
r-dependence of vq for the standard a = 0.17, r s = 20 



kpc profile, normalized so that v = 250 at r 



(the 



latter choice does not effect the shape of r-dependence) 
in figure [SJa). The shape is strongly dependent upon the 
choice of Xj with the larger \ values giving rise to more 
pronounced dependence of vq upon r. 

Our choice of vq (r) differs from one adopted by other- 
authors [63j, i>o(r) ~ r -1 / 4 , that rises as r — » 0. The 
r -1 ' 4 ansatz was inferred from looking at fig. 2 of ref. 
[64J , one of the earlier studies of the effects of baryons on 
the halo properties. (It should be noticed however that 
even in this figure, the logarithmic slope of the z = 
velocity profile does not remain constant as r — > 0, but 
starts to turn downward as in our fig. O^a).) Based on 
the more detailed study [6Gj, we infer that the apparent 
continued rise of vq toward r = is consistent with our 
choice (fT7]l if one also adopts the cuspy profiles found 
by the simulations that include baryons. To illustrate 
this, we plot an example in fig. EJa) (dot-dashed curve) 
using the steepest profile found in the iV-body simulation 
with baryons of ref. [60| . with a — 0.065 and r s = 5.3 
kpc. The upshot is that eq. ([T7| predicts that vo(r) will 
always reach a maximum at some r = r max and thereafter 
fall off as r — > 0, but the value of r max is smaller for 
cuspier halos. In the present work, we will find that 
cuspy halos are not consistent with satisfying the inverse 
Compton constraint, so the i'o(r) ~ r -1 ' 4 would not be 
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the appropriate one for us to use. 

The r-dependence of vq gives rise to r-dependence in 
the velocity-averaged enhancement factor (S) eq. (flO|) 
(la ]. We illustrate for two different values of the cir- 
cular velocity at the solar radius, v c = v a {r & ) = 220 
and 250 km/s, adopting the exponent \ = 1-64 which 
is in the middle of the range for DM simulations includ- 
ing baryons. The enhancement factor is stronger near 
the galactic center than at the Sun, implying that IC 7 
rays are produced more copiously relative to leptons than 
would be the case for a spatially constant boost factor. 
This is illustrated in figure |9jb) where (S(r)) is plotted. 
As a result, the IC bound constrains the models more 
strongly than if one ignored this effect. Like for vo{r), 
the shape of (S(r)) is nearly independent of the value 
of v c , but depends strongly on the choice of exponent %• 
Had we chosen the pure DM value \ — 1-9, we would get 
a much stronger ratio of (S) in the galactic center versus 
solar regions, which would make the IC constraint even 
more difficult to satisfy. 



B. Estimation of lepton and IC 7-ray signals 

Our procedure now is to define the cross section for 
production of e^ (and below, that of IC 7 rays), using 
the partial cross section ((SJ) and enhancement factors that 
have been averaged over the appropriate regions of space. 
The average is done with a weighting factor of r 2 p 2 , ap- 
propriate for annihilations occuring in a shell of radius r 
and thickness dr centered on the sun. Assuming that the 
observed leptons have diffused 110 further than Ar = 1 
kpc [8OJ, the leptonic enhancement factor is 
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(18) 



and the leptonic cross section is a e ± = buB e ±aof, where 
<7o is the reference cross section that would give the cor- 
rect relic density, <r Q v = 3 x 10 -26 cm 3 /s. The resulting 
a e ±v is then compared to the PAMELA and Fermi al- 
lowed regions from [l7|, rescaled by {po,o/p<s) 2 f 2 if nec- 
essary. (We do not apply a correction factor for changing 
a or r s as in the case of the IC constraint, because the 
behavior of p(r) in the vicinity of r© does not change 
significantly as a function of a or r s .) This tells us if 
the predicted leptonic signal is sufficiently close to the 
observed one. 

For the IC 7 ray signal we follow a similar approach, 
but now (S(r)) should be averaged over a radial interval 
that corresponds to the path followed by the relevant 7 
rays. Since the IC signal arrives along a line of sight 
rather than from a fixed volume, the weighting factor is 
p 2 rather than r 2 p 2 . We thus take 
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J^drp 2 (r)(S(r)) 
r°drp 2 (r) 



(19) 



where n = tan(5°)r Q = 0.73 kpc is a cutoff due to the 
fact that the skymap in fig. [3] does not use the inner 5° 
for setting the IC bounds. The predicted 7 ray signal 
corresponds to the cross section 07 = buB^a^f . This 
should be below the bound 
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implied by eqs. (|12I13I15|) . 



C. Search of parameter space 

We have surveyed the predictions of the models for 
a range of DM halo parameters p©, a, r s , v c in search 
of examples that can satisfy all the constraints. There 
are two principal challenges to confront. First, as we 
illustrated in fig. ® the Sommerfeld enhancement tends 
to be too large to fit the lepton signals, especially if the 
mass splitting is nonzero. Second, the IC constraint rules 
out large boost factors unless the DM halo parameters are 
taken to be noncuspy, i.e., large values of a and r s . On 
the other hand, A-body simulations favor smaller values 
of a = 0.16 — 0.17 for galaxies of the size of the Milky 

way mmm. 

The upshot is that a compromise must be made: to 
avoid extremely large values of a, one needs to marginally 
satisfy the IC constraint by making the annihilation rate 
as small as possible, which also pushes the lepton sig- 
nal to the lower boundary of its allowed region. In or- 
der to achieve this reduction in the annihilation rate, we 
are forced to assume that the annihilating DM only con- 
stitutes the fraction 1// of the total mass density, by 
increasing the gauge coupling as in (fj| to suppress the 
relic density. Since the Sommerfeld enhancement rises 
as a function of the mass splitting, there is a trade-off 
between the parameters / and SM. 

We illustrate the best-case scenarios in figure [TU] for 
the triplet DM model with SM = and SM = 1 MeV, 
which shows the model predictions for u e ±v and the 3tr 
confidence level allowed values for PAMELA/Fermi for 
the two cases. In the first case a sufficiently small anni- 
hilation rate is found using the maximum DM fraction 
1// = 0.4, while in the second, it must be reduced to the 
level of 1/7. These examples illustrate the optimal com- 
promise for satisfying the IC constraint; we are able to 
just accommodate it using the large, noncuspy Einasto 
parameter choices a = 0.20, r s = 30 kpc. 6 The lep- 
ton rate is barely consistent with the PAMELA/Fermi 
allowed region for a DM mass of 800 TeV. The relatively 
low circular velocity v c — 220 km/s we have adopted also 
helps to satisfy the IC bound, since as figure [9] implies, 



6 Rcf. 82] also finds a = 0.20 in a recent fit to the data including 
7 rays. 
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FIG. 10: Predicted and allowed values of crv versus M for triplet DM with gauge boson mass /j, — 50, 100, 150 or 200 MeV, and 
halo parameters a = 0.20, r s = 30 kpc, p© = 0.3 GeV/cm 3 . Only predicted values consistent with IC constraint are shown. 
Left: mass splitting SM = 0, DM fraction 1// = 0.4; right SM = 1 MeV, 1// = 1/7. 



the ratio B 1 jB £ 
high v c . 



is somewhat smaller for low v c than for 



To highlight the sensitivity of the IC constraint and 
lepton signal to the cuspiness of the halo profile, we 
present results of our search using the neighboring 
Einasto parameters a — 0.17, 0.25 and the triplet DM 
model, in figure [TTJ The left-hand panels are similar to 
fig. [TTJ1 while the right hand ones plot the ratio of a 1 to 
the IC bound (f2"0l . The points above the dashed line 
in the right-hand panels are ruled out by the IC bound, 
and the figures show how this becomes more severe for 
smaller a. The left panels also show the CMB constraint 
from reionization of the early universe; this is satisfied by 
a much wider margin than is the IC constraint, as already 
anticipated. Figure[TT]illustrates that for the usually pre- 
ferred choice a = 0.17, only compatibility with PAMELA 
can be achieved, but not simultaneously with the Fermi 
lepton anomaly. Only models with M < 700 GeV are 
consistent with the IC constraint, whose strength is an 
increasing function of M. 

The tension we find here is not ameliorated by con- 
sidering the doublet DM model instead of the triplet. 
Given that the gauge couplings a g and cross-section ra- 
tios 6n are rather different between the two models, their 
predictions are suprisingly similar. (This is partly under- 
standable in that we adjust the tree-level cross sections 
of the two models to be the same, but the Sommerfeld 
enhancement introduces extra dependence upon a g , so 
there seems to be an accidental cancellation of the effect 
of 6n by the increase in the Sommerfeld boost.) More- 
over the dependence upon the halo parameters is very 
similar between the two models. In fig. [12] we show the 
predictions and constraints on av versus M for the dou- 
blet model, using the SM = 1 MeV mass splitting, and 
two of the same sets of halo parameters (the less cuspy 
ones) as we considered for the triplet model. The results 
are difficult to distinguish from those of the triplet model. 



VII. DISCUSSION 

We have tried to carefully assess the likelihood that 
some of the best theoretically motivated models of 
annihilating multistatc dark matter can explain the 
PAMELA/Fcrmi/H.E.S.S. lepton excesses, considering 
important details of the models whose implications could 
be missed by a more generic, model-independent inves- 
tigation. As well, we take into account simulations and 
observations of the DM halo properties, which have a 
crucial impact on the viability of the scenario. The 
main challenge to the models is the difficulty of satis- 
fying constraints from Fermi due to the production of in- 
verse Compton (IC) 7 rays produced by the high-energy 
leptons scattering on galactic radiation. 

Among the most important model-dependent features, 
which affect the ease of satisfying the constraints, are 
that each annihilation produces two dark gauge bosons, 
hence four leptons, and these are all e ± as opposed to 
heavier charged particles if the gauge boson masses are 
lighter than ~ 210 MeV. The IC constraints for this case 
are less severe than for models that produce only two lep- 
tons, or ones that produce ^ in addition to e . More- 
over the gauge coupling is bounded from below by the 
requirement that the relic density not exceed what is ob- 
served. We find that the boost factor from Sommerfeld 
enhancement is too large to satisfy the constraints, unless 
the density of the annihilating DM is suppressed by a fac- 
tor 1// < 0.4; hence one would need an additional com- 
ponent of nonleptophilic DM to make up the rest. Using 
recent results for the Sommerfeld enhancement factor of 
multistate DM |38[ , this problem is exacerbated when the 
mass splitting is nonzero; one needs 1// ~ 0.14— 0.2 if 
SM ~ 1 MeV. 

A main result is that the Einasto parameters which de- 
termine the shape of the DM halo must take large values 



0.20 and 



30 kpc in order to barely satisfy the 
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FIG. 11: Left: Predicted and allowed values of av versus M for triplet DM with gauge boson mass fi = 50, 100, 150 or 200 
MeV, SM — 1 MeV, and DM density fraction 1// = 1/7. Only the model predictions that are consistent with the IC constraint 
are shown. Right: ratio of predicted versus maximum allowed IC 7 ray cross section for same models. Top row: a = 0.17, 



25 kpc, and p© = 0.27 GeV/cm J ; bottom row: a = 0.25, 



30 kpc and p Q = 0.35 GeV/cm J 



IC constraint and have marginal consistency with the ob- 
served lepton excesses. In the context of pure dark mat- 
ter simulations, such a large value of a is above the norm 
for a galaxy of the size of the Milky Way. Ref. [81| finds 
a correlation between a and the virial mass of galaxies 
M vir , a = 0.155 + 0.0095 i> 2 , where v = 6 CT it(z)/a(M viT ) 
are functions arising in the Press-Schechter formalism 
83 1 such that 5 cr it = 1-686 at redshift z = 0, and a > 2 
8J] for a galaxy such as ours with M v ; r = 10 12 M Q [70l ]. 
Thus v < 0.843, which implies a < 0.162. However this 
is only true for the mean value, and a fluctuation as large 
as a — 0.20 for a galaxy like ours is possible [83. On the 
other hand, recent results incorporating baryons indicate 
that a is decreased relative to its pure DM value due to 
adiabatic contraction, which exacerbates the conflict. 

We conclude that there is a distinct tension for anni- 
hilating models to consistently account for all the lepton 
excesses, while it is easier to explain that of PAMELA 
alone using smaller DM masses M < 400 GeV. A possible 
loophole that we have not considered here is the sugges- 
tion that leptons originating from subhalos of our main 



halo could provide a substantial fraction of the observed 
lcptonic excess [54j. If this is the case, the production 
of IC 7 rays from the galactic center would be reduced, 
providing another way to weaken the constraint. How 
significant this effect could be is under investigation. 
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Appendix A: Cross section ratios for X1X1 ^ e 

In this appendix we compute the ratio of the cross 
section for ground state DM annihilation into 4e final 
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states relative to the total annihilation cross section in 
the early universe. 



1. Triplet model 

The total cross section for xx ~^ BB annihilation was 
derived for DM in any representation of SU(2) in ap- 
pendix C of ref. [3|| . For the triplet model, the averaged 
squared matrix element can be expressed in the form 



\\M, 



\\-M.\2\ 2 , where \M.%j ] 2 is the matrix element 



for XiXji summed over all possible final states. The fac- 
tor of g comes from averaging over the colors of \j ■ On 
the other hand, the matrix element for XiXi — * B2B2 
was shown to be just half of |A^n| 2 , since the latter in- 
cludes X1X1 ~~ ► B3B3. Therefore the ratio 61 1 for the 
single process X1X1 ~~ * B2B2 is given by 



611 = 



2 l-^ii I 



\\M 



11 



\M 12 \ 



+ 5/2) 7 



(Al) 



using the results of [36[ for the matrix elements. In units 
of g 4 , \Mn\ 2 = 2(|ij 2 + M 2 ) where \t\ 2 and \u\ 2 stand for 
the contributions from the t and u channels, respectively, 



and take the values |t|' 



|iP 



\st\ with \s\ 



- 2. Likewise |^Vf 12 1 2 = : 
-19/4 and \st\ = 4. The 



intefcrence between t and u channels vanishes, \tu\ = 0. 



2. Doublet model 

For the doublet model, the matrix element for XiXj 
B a B{, can be expressed in the form 



cross section for the early universe, we sum over a, b and 
average over i, j. This gives 

1 

4 



(\M\ 2 ) 



S('" 2 



+ \u\ 2 ) + ^\s\ 



12 n , 
T- St 



9 

Tg 



su\) 

(A3) 



where we used the values given above for |i|^ 



This expression agrees with 

_ 33(3+1) 1 



\st\ and also \su\ = \st\ 

the general result (\M\ 2 ) = aj ^> \j(j + 1) - 1/4] for 
the spin-j representation of SU(2) found in [36|], when 
j = 1/2. We need to compare this result with that of the 
exclusive channel X1X1 — ^ B2B2, 



\M{{ 



22|2 



5« 



(A4) 



The ratio of these would give (\)/(jq) = 4/9. However 
unlike the triplet model, there is an additional process 
that can give the B 2 gauge boson, namely X1X1 ~ * B\Bi, 
which has 

1 



\M[{ 



12|2 



16 

1 

Tg 



(|i| 2 + | u | 2 ) + i | s | 2 + -(| s t| + | su |) 

(A5) 



This should be doubled to take into account the X1X1 ~~ ^ 
B2B1 contribution, but we halve it again to account for 
the fact that only one B2 gauge boson is produced, so it 
gives half as many leptons as the B2B2 final state. Hence 
we find that 



&ii = 



_i_ 

16 
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Mt = ~f 



Tuit - 
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tLu + 
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where s,t,u are the spinor parts of the amplitude in the 
respective channel. The factors of 1/2 are for the nor- 
malization of the SU(2) generators. To find the averaged 



Appendix B: Multistate Sommerfeld enhancement 

The Sommerfeld enhancement factor for multicompo- 
nent DM annihilation is different from the that of the 
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single component case. In particular il the gauge cou- 
plings are off-diagonal, Xi-$$X2 an d if there is a mass dif- 



J 



ference 8M between \i an( i X2, ref. 
enhancement factor is given by 



|3g| finds that the 






cosh^-Tr/i;)— cos( We^—e^7r/i/+29- ) 

sh ( (^ + \/ Z 4+^f ) 7r/2;y) soch ( (e„ - y^f+^f) tt/21/) 
cosh((ev+y/-e%+el) ir/iA -cos(20_ ) 



e v < £(5, 



e v > es- 



(Bl) 



r 



with the definitions e v = v/(ca g ) (v being the center-of- 
mass velocity of one of the DM particles), e^ = /i/(Ma), 
and v = |e^(l + y/1 + 4/(e^J"jy)). r« is a radial coor- 

max(e 2 /2,e 2 ), 

where 6,5 



dinate value defined such that V(tm) — 
z g and V{r) = e 



y/25M/M/a g and V(r) = e- d *7r. 6L is a 
function defined as an integral that must be done numer- 
ically (see eq. (4.8) and footnote 3 of [38j). 



Appendix C: Escape velocity 

In general, the escape velocity from a given radius is 
*iven in terms of the gravitational potential 



1 



/OO 773 



(CI) 



where M(r) is the mass enclosed within radius r. If all of 
the matter in the galaxy was dark, we could simply inte- 
grate p(r) to find M(r), but baryons constitute a signifi- 
cant fraction of the matter in the inner part of the galaxy, 
so this would give an underestimate of v esc . Instead, one 
can try to infer the shape of M(r) from rotation curves, 
which measure the circular velocity v c (r) 2 = GM{r)/r. 
We thus have the relation 



dr 



vl(r) 



(C2) 



For flat rotation curves, M(r) increases linearly with 
r, and v^ sc goes like — ln(r). Recent observations of the 



Milky Way 70] indicate that its circular velocity is nearly 
constant to a radius of no = 10 kpc, and falls slowly like 
r p_1 between no and 60 kpc, where p < 1. This corre- 
sponds to M(r) ~ r 2p_1 in the region between 10 and 60 
kpc. Since v c falls from 200 km/s at rio to 170 km/s at 60 
kpc, we infer that p = 0.856. This tells us that v c = 220 
km/s for r < no and v c — (220 km/s)x (r/rio) p_1 in the 
region between no and 60 kpc. 

We still need an estimate of how v c behaves at r > 60 
kpc in order to do the integral (|C2|) . The most radical 
assumption would be that there is no more mass beyond 
this radius, so that v c drops like 1/^/r beyond this point. 
However, the appearance of the rotation curve at the 
highest measured radii does not suggest such a change; 
moreover we have an estimate of the total mass of the 
galaxy through its virial mass M V - 1T , determined to be 
1O 12 M in ref. [70]. Therefore it would seem reasonable 
to assume that M(r) continues to grow with the same 
power law r 2p_1 out to a radius that contains M v i r , after 
which it stops growing. In the present case that radius 
corresponds to r vlr — 175 kpc. We therefore take r v - u as 
our estimate of where M(r) ceases to grow. In the region 
r < no, this gives 
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